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Abstract. Electromagnetic probes promise to be direct messengers of (spectral 
, properties of) hot and dense matter formed in heavy-ion collisions, even at soft 

\Q ' momentum transfers essential for characterizing possible phase transitions. We 

examine how far we have progressed toward this goal by highlighting recent 
I \ developments, and trying to establish connections between lattice QCD, effective 

^ ' hadronic models and phenomenology of dilepton production. 
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O ! 1. Introduction 

■ 

o ■ 

High-energy heavy-ion physics aims at investigating the properties of strongly in- 
teracting matter under extreme conditions, at much higher temperatures (T) and baryon 
'"q ■ densities (qb) than accessible in atomic nuclei. A major goal is to establish the existence 
of new states of matter characterized by the restoration of the spontaneously broken 
chiral symmetry and the deconfinement of color charges, as predicted by Quantum 
Chromodynamics (QCD). Hadronic observables are widely (and successfully) used to 
infer the composition and collective expansion of the highly excited system and thus 
determine its bulk properties, i.e., equation of state (EoS). Direct information on the 
degrees of freedom and spectral properties of the matter are harder to obtain, especially 
at soft scales typical for the phase transition, ~T C , Aqcd, up to the chiral scale A x ~47r/ 7r 
(/ 7r =92 MeV: pion decay constant). Electromagnetic (e.m.) probes (real and virtual 
photons) provide a unique opportunity to study the J p = 1~ response of the strongly 
interacting medium [U EJ [31 0] , defined by the pertinent retarded correlation function, 

lC(g;/^,T) = -ij d'x B(t) (0^(^(0)10) (1) 

finite T, \i B ground state). From e + e~ ^-hadrons one knows that, in the vacuum, 
the low-mass (M< A x ) e.m. spectral function is saturated by light vector mesons, 



imn em ~ 



lmD p + imZ), + ^IimD, 



(2) 



which also illustrates the dominant role of the p. In vector-dominance models, medium 
modifications are thus encoded in the vector-meson propagators, Dy (V = p,u,(p). 
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The study of vector mesons in hadronic matter should be placed in a broader 
context. This includes questions of the type: (1) Do p, to and <fi mesons behave 
alike? (2) Are hadronic interactions sufficient to account for medium effects or are 
"intrinsic" T- and /i^-dependencies of the parameters of the effective Lagrangian to be 
included [SI E]? (3) How do vector-meson spectral functions change across the phase 
diagram (i.e., baryon- vs. meson-induced modifications)? (4) What is the impact of 
medium-modified hadron properties on the EoS and its chemical composition? (5) What 
are the connections between vector spectral functions and (chiral) order parameters 
(quark and gluon condensates, f n , susceptibilities, ...), and how can they be exploited? 
The nonperturbative nature of these questions requires the use of effective (chiral) 
models with careful constraints from both first principles (symmetries, lattice QCD) and 
experiment. Phenomenological applications to electromagnetic observables in heavy-ion 
collisions can then pave the way to deduce the onset and realization of chiral symmetry 
restoration, a key prediction of QCD. 

In this paper, we discuss 3 basic elements (lattice QCD, effective models and pheno- 
menology) to address the above issues. Sec. [2] recalls relations between e.m. spectral 
functions and chiral order parameters via susceptibilities and sum rules. Sec. [3] 
gives a brief survey on vector mesons in hot/dense hadronic matter. Sec. H] assesses 
ingredients for studying dilepton spectra in heavy-ion collisions, with applications to, 
and interpretations of, recent SPS data[jj]. Sec. [S] contains concluding remarks. 



2. Vector Correlators from Lattice QCD and Sum Rules 

The basic quantity which, in principle, contains all the information about an interacting 
system at finite T and \xb (baryon chemical potential) is the thermodynamic potential, 

n(fjL B ,T) = -T\ogZ=-PV, (3) 

defined in terms of the partition function, Z (P: pressure). Its T and /ig derivatives give 
rise to the EoS (entropy, energy-/density), while hadronic correlation functions follow 
from expectation values of suitable currents, cf. Eq. ([!]). The EoS is routinely computed 
in 1QCD, but the formulation of 1QCD in Euclidean (imaginary) time renders it more 
difficult to extract spectral functions in the timelike regime. A valuable interface of bulk 
and microscopic quantities is provided by susceptibilities, which are 2. order derivatives 
of Q that can be related to the spacelike (screening) limit of pertinent correlators, 

Xa~^~n Q ( go = 0,q^0) . (4) 

OfM a 

Fig. [1] shows a recent 1QCD computation of the T-dependence of light quark-number 
susceptibilities in the isoscalar (a=q) and isovector {ct=I) channel, carrying the quantum 
number of the uj and p meson, respectively. With increasing p q (evaluated via a Taylor 
expansion), the isovector channel behaves smoothly while the isoscalar one develops a 

|| Medium effects on vector mesons in eiementary (photo-/electro-/hadro-) production off nuclei have 
been discussed in excellent review talks by Djalali [7] and Metag [5] and will not be reiterated here. 
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pronounced peak. This could be related to the approach to a critical endpoint in the 
phase diagram, possibly induced by the effect of a a-to mixing. 

A powerful tool to relate chiral order parameters to (in-medium) l/-meson spectral 
functions are sum rules (SRs). In QCD-SRs a (subtracted) dispersion relation for a 
correlation function is formulated for spacelike momenta q 2 =—Q 2 <0 [TO] . 

n a (Q 2 ) = n Q (o) + C Q 2 + Q 4 J • (5) 

The spectral function on the right-hand-side is related to quark and gluon condensates 
which arise as nonperturbative coefficients within an operator-product expansion (OPE) 
in 1/Q 2 of U a (Q 2 ) on the left-hand- side. The condensates on the OPE side generally 
decrease in the medium, which, in turn, requires a softening of the spectral function 
under the dispersion integral. The softening can be satisfied by both broadening and 
a downward mass shift [TTJ [T2]. For uo and p mesons the OPE side is very similar 
(governed by 4-quark condensates), while the subtraction constant differs considerably, 
rip(O) = n w (0)/9, implying stronger medium effects (softening) on the p than on the uj. 

For the vector-isovector (p) channel (which dominates the dilepton rate, cf. eq. (TJ])), 
so-called Weinberg and DMO sum rules [TBI [Ti] are of special interest. Using current 
algebra, one can relate moments of the difference between vector and axialvector spectral 
functions to chiral order parameters. In the chiral limit (m 7r =0) one has 

00 , 

/« = - f—s n [Imn y (s) - lmU A {s)} , (6) 

J 7T 


f-2 = f^~F A , f. 1 = f 2 , /„ = (), h = -27ra s (0) (7) 

(tv pion charge radius, F A : radiative pion decay constant, (O): 4-quark condensate). 
Weinberg sum rules remain valid at finite T [15], demonstrating that chiral restoration 
requires degeneracy of the entire spectral functions. Combining 1QCD computations of 
order parameters with effective model calculations of spectral functions thus provides a 
very promising synergy for deducing chiral restoration from experiment [16J. 
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Figure 2. p- meson spectral function in a hadronic many-body approach. Left panel: 
hot hadronic matter [17] at fixed /is =330 MeV (approximately reflecting CERN- 
SPS conditions) corresponding to (g B / Qa,T/MeV)~(0.1, 120), (0.7,150), (2.6,180); right 
panel: hot meson gas (no anti-/baryons) [18] . 



3. Hadronic Spectral Functions 

A commonly used approach to evaluate in-medium V-meson spectral functions is 
hadronic many-body theory pQ. The key quantity figuring into the V- meson propagator, 

D V (M, q; fi B , T) = [M 2 -m\- Z VP - Z VB - E VM ] , (8) 

is the selfenergy, Sy, which can be calculated from effective interactions integrated over 
the momentum distributions of mesons (M=ir,K,p,...) and baryons (B=N,A,A,...) in 
the heat bath (Eyp encodes medium-modified nn, 3tt, KK decays). The interaction 
vertices ought to be compatible with gauge and chiral invariance and constrained by 
resonance decay widths and scattering data in the vacuum. Typical results for the p 
spectral function are shown in Fig. [2j Prominent features are a strong broadening with 
little mass shift and the prevalence of anti-/baryonic effects, especially below the free 
p mass. When extrapolated to the expected phase boundary, the p spectral function 
"melts" and resembles the form and magnitude of the perturbative qq correlator. In 
cold nuclear matter, it is compatible with QCD sum rule analysis [11] and qualitatively 
consistent with the new CLAS data for photoproduction off nuclei at JLab [7] . 

In Ref. [19], the Sv's have been evaluated in terms of empirical scattering 
amplitudes as constructed from resonance dominance at low energies and Regge-type 
behavior at high energies, while the real parts have been calculated using dispersion 
relations. It is satisfying to see that the resulting p spectral function agrees quite 
well with the many-body approach, cf. left panel of Fig. [3] (the largest discrepancy 
emerges toward the 2tc threshold; also, some care has to be taken in identifying nucleon 
and baryon densities at finite T). The u spectral function (right panel of Fig. [3]) also 
broadens considerably (with little mass shift), but less than the p, which is in line with 
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Figure 3. y-meson spectral functions in a hot tt-N gas based on empirical scattering 
amplitudes constrained by resonances and Regge behavior [19] . Left panel: p meson, 
compared to hadronic many-body calculations of Fig. [2] Right panel: uj meson. 



the expectation from QCD sum rules (recall the discussion after eq. 

Within the Hidden Local Symmetry approach, it has been suggested to identify 
the chiral partner of the (longitudinal) p with the 7r (rather than a\) [5j. When 
applied within a finite-T loop expansion, and after matching the hadronic axial-/vector 
correlators to pQCD in the spacelike regime, one finds a dropping p mass as well as 
a violation of vector dominance. The finite-T e.m. formfactor, which enters into the 
dilepton production rate, clearly exhibits the downward moving p peak [20]. It would 
be interesting to see whether this feature persists when including the effects of baryons. 



4. Dilepton Phenomenology in Heavy-Ion Collisions 

The dilepton emission rate is directly proportional to the e.m. spectral function [2"T] . 

(/A " f B (qo;T)lmU em (M,q;p B ,T) . (9) 



d 4 xd 4 q tt 3 M 2 

Assuming the medium in a heavy-ion collision to kinetically equilibrate, thermal dilepton 
spectra can be calculated by convolution over a realistic space-time model. Obviously, 
the latter must be consistent with hadronic chemistry (evolution of pb and T) and px 
spectra (collective expansion). However, depending on centrality and invariant mass or 
of the lepton pair, nonthermal effects/sources, such as Drell-Yan annihilation and 
final-state decays (including correlated open-charm), may not be neglected. 

In the upper panels of Fig. H] thermal fireball calculations employing V^-meson 
spectral functions from hadronic many-body theory [22J are compared to recent NA60 
data (where improved experimental accuracy enabled a subtraction of final-state hadron 
decay contributions, except for p and charm decays) [23J. The low- mass region is 
well described by the in-medium spectral functions (upper left panel), supporting the 
predicted p line shape (QGP emission is rather small). When switching off baryonic 
effects (or all medium effects), the spectra are incompatible with the data (upper right 
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Figure 4. NA60 dimuon invariant-mass spectra in central 158 AGeV In-ln 
collisions [23], compared to thermal fireball emission based on hadronic many-body 
theory and qq annihilation in a QGP 22J (upper panel; left: full results, right: 
without baryon effects and when using the free e.m. spectral function). Lower panels: 
thermal spectra using e.m. rates from a chiral virial expansion within a hydrodynamic 
evolution 26J (left) and a thermal fireball [22] (right). 



panel). In fact, the space-time integrated dilepton line shapes for different scenarios 
reflect the changes of the p spectral function in Fig. [2], confirming good sensitivity to 
medium effects. At masses above 1 GeV, the free e.m. correlator mostly couples to 
4-pion states. Inclusion of this contribution in the fireball emission leads to reasonable 
agreement with the observed enhancement at M>1 GeV (effects of leading-T vector- 
axialvector mixing [24 J increases the 4tt yield slightly). Note that the relative strength 
of the various thermal sources is fixed; the absolute yield is sensitive to the fireball 
lifetime, which in Fig. H] amounts to about 7 fm/c (the initial/final temperatures are 
7^0^195, 120 MeV). The lower panels of Fig. H] display results of using the in-medium 
e.m. spectral function calculated in a 3-flavor chiral reduction formalism [25] , which 
couples chiral Ward identities with a low-density expansion using empirical p-fj-N 
and p-7r scattering amplitudes. The low- and intermediate-mass enhancement is well 
described, but the lack of p-broadening (inherent to the low-density expansion) leads to 
an overestimate around the free p mass. The spectra obtained from the hydrodynamic 
evolution used in Ref. [26J (left panel) agree well with the fireball approach of Ref. [22] 
(right panel). This also holds for the (smallness of the) perturbative QGP contribution. 
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Figure 5. Dimuon qr spectra in semicentral 158 AGeV In-In collisions [28] for 
mass bins 0.6 GcV<M<0.9 GeV (left) and 1.0 GeV<M<1.4 GeV (right), compared 
to thermal fireball calculations with in-medium p, lo, <f> spectral functions, An 
annihilation, QGP emission, thermal-freezeout and primordial p's, as well as Drell- 
Yan annihilation 1311. 



On the contrary, in Ref. [27] QGP radiation is the dominant source above M ~ 0.9 GeV. 

Additional information on the nature of the excess radiation in heavy-ion collisions 
can be obtained from pair-g^ spectra, especially when differential in mass, see Fig. [5j At 
low qx the spectra are dominated by thermal radiation (either 2ir or 4tt annihilation), 
but the composition becomes more involved at higher qr- In the p-mass region, p decays 
at freezeout (subject to a large flow-induced blueshift, with in-medium spectral shape) 
as well as primordial p's from minijets escaping the medium (including a Cronin effect, 
with vacuum spectral shape) become significant. Note that the emission kinematics of 
these sources actually differs from thermal radiation by a time dilation factor r )=q§fM. 
In the p-mass bin, the theoretical qx spectra appear somewhat too soft (as more clearly 
revealed in a local slope analysis [28]), cf. also Ref. [27]. An interesting possibility 
is that viscosity effects in the later hadronic stage render hadron spectra harder [29J. 
Transverse-momentum and centrality dependencies certainly deserve further study. 

The CERES/NA45 collaboration has recently published dielectron excess spectra 
for central Pb(158 AGeV)+Au [30]. Albeit inferior in statistics to the NA60 dimuons, 
the small electron mass allows access to the mass region below the 2tc threshold, where 
theory predicts a large enhancement (the data might show first indications of this). 
In addition, in larger collision systems, the relative size of nonthermal contributions is 
suppressed. Hadronic many-body calculations describe the CERES spectra well. 

5. Conclusions 

E.m. probes are progressively improving our understanding of microscopic properties 
of hot and dense matter. In this paper we have focused on soft emission, which 
predominantly relates to the physics of in-medium vector mesons. Hadronic models, 
when carefully constrained, predict a strong in-medium broadening ("melting") of the 
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p resonance which has withstood a substantial increase in accuracy in recent dilepton 
measurements at the SPS. This also corroborates the notion of (anti-) baryons as the 
major agents of medium effects, which is expected to hold even at RHIC and LHC, 
raising intriguing questions on the nature of chiral restoration. The determination of 
the in-medium axialvector spectral function remains a pressing issue. It can and will 
be addressed by constructing realistic chiral models in a hot and dense medium, which 
will serve as a bridge between order parameters (as computed, e.g., in lattice QCD) and 
experimental data, by utilizing chiral and QCD sum rules. 
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